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ABSTRACT 

We have analyzed solar-like oscillations in ^1700 stars observed by the Kepler Mission, spanning 
from the main-sequence to the red clump. Using evolutionary models, we test asteroseismic scaling 
relations for the frequency of maximum power (v max ), the large frequency separation (Az/) and os- 
cillation amplitudes. We show that the difference of the A^-^ max relation for unevolved and evolved 
stars can be explained by different distributions in effective temperature and stellar mass, in agree- 
ment with what is expected from scaling relations. For oscillation amplitudes, we show that neither 
(L/M) s scaling nor the revised scaling relation by Kjeldsen & Bedding (2011) is accurate for red-giant 
stars, and demonstrate that a revised scaling relation with a separate luminosity-mass dependence 
can be used to calculate amplitudes from the main-sequence to red-giants to a precision of ~25%. 
The residuals show an offset particularly for unevolved stars, suggesting that an additional physical 
dependency is necessary to fully reproduce the observed amplitudes. We investigate correlations be- 
tween amplitudes and stellar activity, and find evidence that the effect of amplitude suppression is 
most pronounced for subgiant stars. Finally, we test the location of the cool edge of the instability 
strip in the Hertzsprung-Russell diagram using solar-like oscillations and find the detections in the 
hottest stars compatible with a domain of hybrid stochastically excited and opacity driven pulsation. 
Subject headings: stars: oscillations — stars: late-type — techniques: photometric 
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1. INTRODUCTION 

Empirical relations connecting observable quantities 
with physical parameters of stars are of fundamental im- 
portance for many fields of stellar astrophysics, with clas- 
sical examples including the colour-temperature scale for 
late- t ype stars (see, e.g.. lFlowerlll996t ICasagrande et al] 
I201CO an d the period-l uminosity relation for Cepheid 
variables (|Leavittl H908). The calibration of these rela- 
tions relies either on the direct measurement of stellar 
properties (e.g. through trigonometric parallaxes and 
interferometry) or the comparison of observations with 
stellar models. In both cases measurements spanning an 
extended range in parameter space (such as temperature, 
luminosity and metallicity) are required. 

One of the most powerful methods to constrain funda- 
ment al_j£op_ertfo^_ci_field stars is asteroseismology (see, 
e.g., [ Brown fc Gillilandl 119941 : iChristensen-Dalsgaardl 
120041: lAerts et al.1 12010 ). The connection of global 
asteroseismic observables t o fundamental stella r prop- 
erties was calibrated by iKieldsen fc Bedding] (|1995l 
hereafter KB95) by introducing relations that scale 
stellar properties from the observed values of the 
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Sun. The scaling relations have been used extensively 
both for the direc t inference of stellar properties (e.g., 
Steilo et al.l l2009bl: iKallinger et al.ll2010at IMosser et all 



20ll" 



Hekker et al.1 



Kallinger et al l l2010bt 
l2011bl lal: 



Chapli n et all [2011a; 
iSilva Aguirre et al.1 120 111 ) 
as well as for calculating pulsation a mplitudes of 
main -se quence and su b giant st ars (e.g.. IMichel et all 
20081: iBonanno et all 120081: iMathur et all l2010bl: 



Huber et al. 2011: iChaplin et all 1201 id) and re d -giant 
stars (lEdmonds fc Gillilandl 119961: iGillilandl I20F 
Stello fc Gillilandll2009t IMosser et al 1120101 iHuber et al 



20101: IBaudin et all 1201 li IStello et al.l l27)Tlj) . The three 



global observables discussed in this paper are the fre- 
quency of maximum power (^ max ), the large frequency 
separation {Av) and the mean oscillation amplitude per 
radial mode. 

Tests of the scaling relations for these observables over 
a wide range of parameter space has so far been ham- 
pered by the relative sparsity of detections of solar- 
like oscillations. The success of the Kepler mission has 
changed this picture, with the detection of oscillations in 
thousands of stars covering a large part of the low- mass 
region in the H-R diagram (|Gilliland et al.l l2010bl ). In 
this paper, we aim to use the large number of detections 
by Kepler to test scaling relations for stars ranging from 
the main-sequence to He-core burning red-giants. 

2. DATA ANALYSIS AND MODELS 

The Kepler space telescope was launched in March 
2009 with the primary goal of finding Earth-like planets 
orbiting solar-like stars through the detection of photo- 
metric transits. Kepler monitors the brightness of stars 
at two sampling rates, either in 29.4 min (long cadence) 
or in 58.8 sec (short cade nce) intervals (Gillilan d et all 
l2010al Denkins et al.ll2010h . For asteroseismic studies of 
solar-like oscillators, the former are primarily used to 
study oscillations in red-giant stars, and short cadence 
data are used to measure the more rapid oscillations in 
main-sequence and subgiant stars. While red giants ob- 
served in long cadence have been continuously monitored 
since the launch of the mission, the number of short ca- 
dence slots is restricted due to bandwidth limitations. 
Short cadence slots have therefore so far been primarily 
used to survey a large number of stars for a period of one 
month each. 

Kepler observations are subdivided into quarters, 
starting with the initial commissioning run (10 d, 
Q0), followed by a short first quarter (34 d, Ql) and 
subsequent full quarters of 90 d length. Our studies 
are based on Kepler data spanning from Q0 to Q6 for 
long cadence data and Q0 to Q4 for short cadence data. 
Our final sample contains 1686 stars, of which 1144 
have been observed at long cadence, mostly since the 
launch of the mission (~ 500 days), and 542 at short 
cadence for a typical length of one month. We have used 
Kepler raw data, which were re duced in the manner 
described by iGarcia et all (|201lD and analyzed using: 
sever al automated ana l ysis m e thods (IBonanno et all 
20081: ICampante et all 120101: IHekker et all^ OlOb: 
Huber et al.l 120091: IKallinger et al.l l2010ab iKaroff et all 
^201 Oat 



IMathur et all 



Mosser fc Appourchai 



20K 
200' 



„ IMosser et all I2011bt IVerner fc Roxburgh! 120111) . 
We refer the read er to IHekker et all (|2011d) and 



IVerner et al.l (|2011aD for an extensive comparison of the 



results provided by these methods. Unless otherwise 
mentioned, all results p r esente d here are based on the 
method by IHuber et al.l (|2009l) . We have only retained 
results for stars with at least one matching pipeline 
result within 10% and 5% of the determined f max and 
Av value, respectively. The same outlier rejection 
procedure was repeated for all methods that returned 
results for both long cadence and short cadence data, 
and these datasets have then been used to validate all 
results and conclusions reported in this paper. Note that 
all amplitudes shown in this paper have been normalized 
to recover the full sine-amplitude of an injected signal 
(commonly referred to as peak-scaling). All amplitudes 
sh own have been calcula ted using the me thod described 
bv iKieldsen et all (|2008j ) with c = 3.04 (jBedding et all 
to convert to amplitude per radial mode. 

Uncertainties on i/ max , Av and amplitudes reported 
in this paper were estimated using Monte-Carlo simu- 
lations by generating synthetic power spectra following 
a x 2 distribution with two degrees of freedom with ex- 
pected values corresponding to the observe d power den- 
sity le vels. For each star, the method by IHuber et all 
(2009) was repeated on each synthetic dataset and the 
standard deviation of the distribution after 500 iterations 
was taken as an estimate of the uncertainty. The typi- 
cal relative uncertainties obtained using this method for 
long cadence and short cadence data are 3% and 4% for 
f-Wx, 1% and 3% for Av, and 7% and 11% for the ampli- 
tude. These estimate s agr ee with the results ob tained by 
IHekker et afl (|2011cT) and IVerner et all (|2011a| ). For our 
analysis, only stars with uncertainties lower than 20% in 
f-Wx, 10% in Av and 50% in amplitude were retained. 
Note that oscillation amplitudes for all stars observed in 
long-cadence with j/ max > 200/iHz have been removed 
from our analysis due to the difficulty of estimating the 
noise level close to the Nyquist frequency of 283/iHz. 

The asteroseismic relations discussed in this paper rely 
on scaling from observed values of the Sun. To ensure 
that these reference values are consistent with our analy- 
sis method, we used 111 30- day subsets of data c ollected 
by the VIRGO instrument ()Fr5hlich et al.l[T997l ) aboard 
the SOHO spacecraft spanning from 1996 to 2005 and 
analyzed them in the same way as the Kepler data. This 
yielded solar reference values of t-Wx,© = 3090 ± 30^Hz 
and Avq — 135.1 ± 0.1/iHz, which are consistent with 
previously quoted values in the literature. We deter- 
mined the oscillation amplitude to be A Q = 4.4±0.3ppm 
for the VIRGO green channel (A = 500 nm) which, us- 
ing the approximation by KB95, translates into a solar 
bolometric amplitude of Aq^ \ = 3.5 ± 0.2 ppm. This 
is in excellent agreement with the solar refere nce value 
of 3.6 ppm established bv IMichel et all ([20091 ). and we 
have adopted Aqtq Q \ = 3.6 ppm for the remainder of this 
paper. 

In the absence of directly determined fundamental 
properties for most of the Kepler stars, a test of as- 
teroseismic scaling relations relies on the comparison 
of observations with models. In o ur study, we use 
canon ical BaSTI evolutionary models ([Pietrinferni et all 
120041) with a solar-scaled distribution of heavy elements 
(jGrevesse fc Noels! [19931 ) . Mass loss in BaST I models 
is ch aracterized according to the Reimers law ([Reimersl 
119751 ). and we used models with the mass loss parame- 
ter set to the commonly used value 77 = 0.4 (see, e.g., 
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Fig. 1. — Upper panel: Av versus f ma x for the entire sample of 
Kepler stars. Red triangles show stars observed in long cadence, 
while black diamonds are stars observed in short cadence. Lower 
panel: Same as upper panel, but with the luminosity dependence 
removed by raising z; m ax to the power of 0.75. Green lines show 
power law fits to the Ai/-y max relation for two different intervals 
of t/max (se e text). The blue dashed line shows the relation by 
Stcllo ct al. (2009a). 



iFusi-Pecci fc Renzinil [l976h Note that mass loss for 
these models is only significant in the red giant phase 
of stellar evolution. 

3. SCALING RELATIONS FOR z/ max AND Av 

iBrown et al.l (1991) first argued that the frequency of 
maximum power (i/ mra ) for Sun-like stars should scale 
with the acoustic cut-off frequency. KB95 used this as- 
sumption to relate v max to stellar properties as follows: 



M/M Q (T cS /T cS , 
L/L Q 



\3.5 



" ^max,0 



(i) 



This scaling relation has since b een found to work 
well both observational l y (see, e.g.. IBedding fc Kieldsenl 
[20031: IStello et al.l l2008HBedding]|201lD as well as theo- 
retica lly (see, e.g.. IChaplin et aL 2008; IBelkacem et al.l 

unni). 

The mean large frequency separation (A^) between 
modes of consecutive radial overtone and equal spher- 
ical degree is directly related to the sound travel time 
across the stellar diamet er, and is th erefore sensitive to 
the mean stellar density (|Ulrichlll986f ). This is expressed 
in the following scaling relation: 



Av 



(M/M Q ) - 5 (T eB /T e g, e y 



(L/L e ) 



0.75 



Av G 



(2) 



Note that in our analysis Av is measured as the mean 
spacing of all detectable modes around the value of f max 
in the power spectrum. 

It has been well established for both main-sequence 
and red-giant stars that ^ max and Av follow a power 



TABLE 1 

Coefficients of the Ay-y ma x relation. 



Method 


a,/3 


a,/3 


# of stars 




fmax < 300/lHz 


fmax > 300/^Hz 




A2Z 


0.259(3),0.765(2) 


0.25(1), 0.779(7) 


919,257 


COR 


0.267(3),0.761(2) 


0.23(1),0.789(6) 


1150,415 


OCT 


0.263(3),0.763(2) 


0.20(1),0.811(7) 


1082,281 


SYD 


0.267(2),0.760(2) 


0.22(1), 0.797(5) 


1228,458 



A2Z - IMathur et alJ H2010al) . CO R - IMosser k, Appourchauxl 
(2009); Mos ser et al.l Il2011bl) . OCT - IHekker et al.1 Il2010bl) . SYD - 
IHuber et al.l l|2009l ). 



law relation (IStello et al.1 I2009at IHekker et al.1 120091: 
IMosser et aD2010HHekker et al.H2011bUal) : 



Au = a(f max //iHz) 



(3) 



Figure 1 shows this relation for the entire Kepler sample 
in our analysis. Although the relation a ppears to be 
consta nt ov er several orders of magnitude, Moss er~et al.1 
(|2010fl and IHuber et all (|2010l ) noted that the slope is 
different for red-giant and main-sequence stars. This can 
be illustrated more clearly by removing the luminosity 
dependence by raising j/ max to the power of 0.75, yielding 



(^max/^Hz) 

Av/fiUz 



0.75 



X 



M 



T e S,E> 



-0.375 



(4) 



The lower panel of Figure 1 displays this ratio as func- 
tion of z^ max . The distribution shows a prominent di- 
agonal structure around ^ max ~ 20 — 50/iHz, which we 
identify as the red c lump, comprising He-core burning 
red gi ant stars (see IBedding et all 12011: M osser et al.1 
l2011allH ). It is evident that the power law becomes 
steeper as ^ max increases. To measure this effect, we fit- 
ted Equation © to the sample in two groups subdivided 
at f max = 300/xHz, which roughly marks the transition 
from low-luminosity red giants to subgiants. The best- 
fitting power laws for the SYD pipeline are shown as 
solid green lines in the lower panel of Figure [TJ and the 
results for each pipeline with values for both long- and 
short cadence data are listed in Table 1. The difference 
in the coefficients between red giant and main-sequence 
stars is significant, and should be noted when using the 
relation for determining Av from f max . For i/ max close 
to the solar value, for example, the use of a power-law 
relation calibrated to red-giant stars would lead to an 
underestimation in Au by ~ 10%. For comparison, the 
blue dashed line in the lower panel of Fi gure [T] shows the 
relation derived by IStello et al.l (|2009a[) using a sample 
including both main-sequence and red-giant stars. 

How do our observed values of v max and Av com- 
pare with evolutionary models? The upper panel 
of Figure 2 compares the observed distributions with 
solar- metallicity models with masses l.OAf©, 1.3Af Q 
and 2.0M Q (solid lines), which roughly correspond to 
the lower bound, median an d upper bound of the 
mas s distribution derived by Lkallingc r et al.l (|2010af ) 
and IChaplin et all (|2011a|) . The model values for ^ max 
and Av have been calculated using Equations (HJ and 
([2]). Additionally, we have color-coded the mass of 
each star calculated using Equations ([T]) and ([2]) with 
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100 1000 

Fig. 2. — Upper panel: Observed values of u^l^/ Av versus i^max 
(symbols) compared to solar-metallicity CY = 0.273, Z = 0.0198) 
models with lMg, 1.3Mq and 2.0Mq (solid lines). Asteroseis- 
mic masses are color-coded as indicated in the plot. Typical error 
bars for different ranges of f max are indicated near the bottom of 
the plot. Lower panel: Same as top panel but with the effective 
temperature dependence removed (see Equation ) and omitting 
model tracks. Symbol types and colors are the same as in Figure 
[T] The green solid lines show linear fits to the linear-log plot for 
the same intervals of i/ ma x as in Figure [T] 



the observed values of v max and Ais adopting the ef- 
fectiv e temperature liste d in the Kepler Input Catalog 
(KIC. lBrown et al.ll^Oill) 1 . The observed and theoretical 
masses in this plot essentially correspond to a compari- 
son of the so-called direct method and grid-bas ed method 
of es timating asteroseismic masses (see, e.g., IGai et all 
I2011D . assumin g solar metallici t y. 

As n oted by iKallinger et al.l |2010b[ ) and iHuber et all 
(2010), the spread in stellar mass is pronounced on the 
red giant branch, while for less-evolved stars the spread 
in observations is weaker and the models almost over- 
lap. The overall agreement between the models and the 
data is very good, and we do not observe any signifi- 
cant offset of the models with respect to the observa- 
tions. It is also remarkable how well the models of dif- 
ferent masses track the He-core-burning red clump. Al- 
though the spread of data points about the models for 
main-sequence stars can be explained by measurement 
uncertain ties, we note that ste l lar po pulation models pre- 
sented by iSilva Aguirre et al.l (|201lP ) yield evidence that 
the Kepler sample is on average metal-poor, which can 



have a some impact on the AiA i w 
sequence stars (see 
ISilva Aguirre et all 



relation of main- 

Ii 1 1 1 1 fl< 1 1 ■ 
12010( 1 . As noted by 

20111 ). however, the measurement 



1 Note that for 39 stars in our list, no effective temperatures were 
available in the KIC, and we omitted those stars for the remainder 
of our analysis. 



uncertainties are currently too large to test a metallicity 
influence on the scaling relations. 

What causes the different Ai/-i/ max power-law relations 
for evolved and unevolved stars? The scaling relation in 
Equation (j4j suggests that the change in slope for the 
unevolved stars must be partially due to a variation in 
effective temperature. To test this hypothesis, we again 
used the effective temperatures listed in the Kepler Input 
Catalog to model this dependency using Equation ©. 
The result is shown in the lower panel of Figure 2. As 
expected, correcting for the higher average effective tem- 
peratures of main-sequence stars compared to red giants 
removes the gradient in the distribution. Subdividing the 
sample again at f max = 300/iHz, we find a positive slope 
for the red giant sample and a negative slope for main- 
sequence stars (solid green lines). These variations are 
qualitatively in agreement with different mass distribu- 
tions in the sample: while for red giants i/ ma x is correlated 
to stel lar mass for He-core burning stars (jMosser et al.l 
l2011allcT ). main-sequence stars wit h higher ^ max are gen - 
erally low-mass stars (see, e.g., iChaplin et al.l l2011al) . 
Without mass estimates independent of asteroseismic 
scaling relations, however, it is not possible to make fur- 
ther quantitative conclusions about the distribution. 

4. AMPLITUDES 

4.1. The L/M Scaling Relation 

KB95 suggested that mod el pre dictions by 
IChristense n-Dalsga ard fc Frandsenl ()1983D implied a 
scaling for velocity amplitudes of 



A vel oc 



M 



(5) 



with 8=1. They further argued that the oscillation 
amplitude A\ observed in photometry at a wavelength A 
is related to the velocity amplitude: 



A x oc 



(6) 



The exponent s has since been revised i n the range of 
roughly s = 0.7 — 1.3 both theoretically ([Houdek et al.l 
Il999al iHoudekl l200l iSamadi et ail l2007h. as well as 
observationally using red - giant st ars (iGillilandl 1 20081 : 
Dziembowski fc Soszvriskil (20101 : iMosser et aiT bold 
Stello et al.l l20l6T T main- sequence stars (jVerner et al.l 
201 laj) and an ensemble o f main-sequence and red-giant 
stars (jBaudin et al.l l20Tll) . The value for r has so far 
been chosen to be either r = 1.5 (assuming adiabatic os- 
cillations) or r = 2.0 (the best-fitting coefficient found 
by KB95 for classical pulsators). 

Since our sample spans a large range in effective tem- 
perature, we must account for the spectral response of 
the Kepler bandpass to compare predictions made with 
Equation to our observations. To do so, we used the 
expression for the bolomctric amplitude given by KB95 
and converted these values to amplitudes observed in the 
Kepler bandpass, as follows: 



.4 



Kp 



1 



T^ s c K (T ei 



(7) 



with ck being the bolometric correction factor as a 
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functi on of effective temperature, given by iBallot et al.l 
(|20T1 : 



ck(T ( 



off J 



cif 



5934K 



(8) 



The reason for choosing to correct the model amplitudes 
rather than the observed values is that effective tempera- 
tures for most of the stars in our sample are rather uncer- 
tain, and we hence prefer not to perform the correction 
on the observed amplitudes. 

Figure Ejja) shows the observed amplitudes for the 
full Kepler sample as a function of v max . The spread 
is much larger than t he typical measurem ent uncertain- 
ties. As fir st noted bv | Huber et al.l (I2010f) and la t er con- 
firm ed bv IMosser et i all (|2011aD , iMosser et~aT1 (|2011cD 
and iStello et all (|20lilK this spread in the amplitude- 
fmax relation for red giants is r elate d to a spread in 
mass. To demonstrate this, Figure l4~lT b) shows the same 
plot but color-coded by stellar masses, as calculated in 
the previous section. We observe that, particularly for 
low-luminosity red giants, the higher-mass stars show 
lower amplitudes than lower-mass stars for a given ^ max - 
For unevolved stars we can tentatively identify the same 
trend, although the separation is less clear. 

The observation of a mass dependence for a given ^ max 
indicates that Equation (|7|) should be revised to include 
an additional mass dependence. To demonstrate this, 
solid lines in Figure I4.1f b) show model tracks with dif- 
ferent masses scaled using typical values of r = 2 and 
s = 0.8 in Equation 0. The (L/M) s scaling clearly 
fails to reproduce the observed spread on the red-giant 
branch, but predicts a strong mass dependence for un- 
evolved stars, contrary to what is observed. 

An obvious way to account for an additional mass de- 
pendence is to rearrange Equation ([7]) as follows: 



.4 



Kp 



(9) 



Note that such a fo r mulat ion has been introduced by 
IKieldsen fc Bedding! d20H ( hereafter KB11) and has 
also been used bv lStello et al.l (|201lD in the Kepler study 
of cluster red giants. To evaluate the coefficients that 
best reproduce our observations, we again used models 
with masses l.OM , 1.3M and 2.0M Q and calculated 
their expected amplitudes according to Equation @ for 
a given set of s, t and r. The agreement of models and 
observations was evaluated as follows: For each observed 
amplitude, we interpolated each model track to obtain 
the model amplitude at the observed f max value. We 
then calculated the minimum squared deviation of the 
observed amplitude to the three model amplitudes, nor- 
malized by the measurement uncertainty. The parame- 
ters s, t and r were then optimized using a least-squares 
fit. Note that for simplicity we have only used model 
values below the tip of the red-giant branch, since an 
inclusion of later evolutionary stages would require the 
identification of all red-clump stars in our sample which 
is beyond the scope of this paper. 

First tests quickly showed that the strong correla- 
tion between effective temperature and luminosity made 
an indepe ndent determinatio n of s, t and r impossible. 
Following Ste llo et all (|20T1 . we have therefore fixed 




Fig. 3. — (a) Oscillation amplitude versus v ma x for the entire 
Kepler sample. Symbol types and colors are the same as in Figure 
1. The position of the Sun is also marked, (b) Same as panel (a) 
but with asteroseismic masses color-coded. Error bars have been 
omitted for clarity. Solid lines show the 1.0, 1.3 and 2.0 Mq solar- 
metallicity (Y = 0.273, Z = 0.0198) models scaled using Equation 
(0 with typical values of r = 2 and s = 0.8. (c) Same as panel (b) 
but using Equation (|9j with coefficients determined by fitting the 
model tracks to observations, (d) Same as panel (c) but with best 
fitting coefficients determined by comparing observed amplitudes 
with calculated amplitudes using asteroseismic masses and radii 
(see text for details). 
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Fig. 4. — Observed versus calculated amplitudes for the best- 
fitting coefficients in Equation (0, using radii and masses deter- 
mined using the scaling relations for ^ ma x and Au. The dashed 
blue line shows the 1:1 relation. The position of the Sun is also 
shown. 



the value to r = 2. The best-htting parameters are 
s = 0.886 ± 0.002 and t = 1.89 ± 0.01, with uncertainties 
estimated by repeating the fitting procedure 1000 times 
using amplitudes drawn from a random distribution with 
a scatter corresponding to the 1 a measurement uncer- 
tainties (scaled so that reduced \ 2 — 1 f° r the original 
data). The parameter uncertainties were then estimated 
by calculating the standard deviation of the resulting dis- 
tribution for each coefficient. 

The scaled model tracks in Figure I4.1f c) using these 
coefficients reproduce a mass spread in amplitude for 
red giants, with higher-mass stars showin g lower am- 
plitud es for a given t- Wy, as observed by jHuber et al.l 
d201Ch. IMosser et al l (1201 lali. IStello et all (|2011h and 
iMosser et al.l ( 201 Id ). For unevolved stars, the model 
amplitudes now show less dispersion but appear to sys- 
tematically underestimate amplitudes for both subgiant 
and main-sequence stars. 

The method presented above assumes that the excess 
spread of amplitudes at a given j/ max is entirely due to a 
spread in stellar mass. A more sophisticated approach is 
to account for the actual mass distribution by directly us- 
ing masses and radii of the sample estimated using z^ max 
and Av in Equations ([TJ and (j2J|. We again used the 
effective temperatures from KIC and evaluated Equa- 
tion © for each value of s and t directly using these 
asteroseismic masses and radii. The corresponding best- 
fitting parameters in this case were s — 0.838 ±0.002 and 
t = 1.32 ± 0.02, and the scaled model tracks with these 
parameters are shown in Figure l4~TT d). The spread of 
the model tracks on the red giant branch is now consid- 
erably reduced, better matching the observations, and 
the overall fit for less evolved stars is also improved. 

Figure [4] compares the observed and calculated ampli- 
tudes using s — 0.838 and t = 1.32, showing good agree- 
ment along the 1:1 line over the entire range of ^ max - 
The residuals between observed and calculated ampli- 
tudes over the full range of f max show a standard devi- 
ation of about 25%, which is consistent with the typical 
uncertainties in the adopted stellar properties (estimated 
from propagating the uncertainties in i/ max and Av and 



assuming an uncertainty of 200 K in T e g) and the mea- 
surement uncertainties in the observed amplitudes. We 
do see an overall bias of about 9%, with calculated am- 
plitudes being systematically underestimated when scal- 
ing from the solar value. This bias is stronger for un- 
evolved stars with ^ max > 300/iHz (15%) than for evolved 
stars f max < 300/iHz (7%). An explanation for this bias 
might be additional physical differences between the Sun 
and our Kepler sample, which influence oscillation am- 
plitudes but have not yet been taken into account. 

Our best-fitting values for s and t using a steroseismic 
masse s and radii are significantly different to IStello et al.l 
poll , who found s = 0.90±0.02 and t = 1.7±0.1 using 
independently determined properties of cluster red giants 
in the Kepler field. This difference is presumably due to 
the fact that our sample includes a much wider range of 
evolutionary states. Indeed, repeating our analysis using 
only red giants in a ^ max range similar to that used by 
IStello et alj (|2011l ) yields coefficients which are in bet- 
ter agreement. Remaining differences could potentially 
be due to metallicity effects which have been suggested 
to have a significant influence on oscillation amplitudes 
(|Samadi et al.l[2010l ) or systematic errors when estimat- 
ing stellar properties through the direct method of using 
f-max and Av. Nevertheless, all methods confirm that 
scaling relations with a separate mass and luminosity de- 
pendence better reproduce the observed amplitudes from 
the main sequence to red giants. 

4.2. The KB 11 Scaling Relation 

Kicl dsen fc Beddinel (|201lD recently argued that am- 
plitudes of solar-like oscillations should scale in propor- 
tion to fluctuations due to granulation. They proposed 
a revised scaling relation for velocity amplitudes: 



-4 V 



Lt, 



0.5 



J oc 



(10) 



where r osc is the mode lifetime. 

Using the same arguments as in the previous section, 
the KB 11 relation predicts photometric amplitudes in 
the Kepler bandpass as follows: 



.4 



Kp 



r -7-0.5 



Ml .5 Tc l ff 25 + , Ck(Teff) 



(11) 



They also suggested the following relation for the gran- 
ulation power in intensity measured at ^ max : 



M 3 T 5 « 5 

err 



(12) 



First results on granulation prope rties of Kepler red- 
giant stars by [Mathur et al.l (1201 ID have shown pr omis- 
ing agreement. More recently. IMosser "eTall ll20llch con- 
firmed the proportionality between mode amplitudes and 
the granulation power at f max for a large sample of red 
giants, but found that the absolute values calculated us- 
ing Equations (fTTj) and (IT21 were offset from the observed 
values. Here we extend this comparison to include also 
main-sequence stars. 

To test the relation, we used the solar granulation 
power at ^ max determined by our analysis of VIRGO 
data, which yielded a mean value of Pint,© = 0.19 ± 



Testing Asteroseismic Scaling Relations 



7 
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Fig. 5. — Upper panel: Obs erve d versus calculated granulation 
power at f ma x using Equation i l 1 21 ). Lo wer panel: Observed versus 
calculated amplitudes using Equation mil l and assuming a solar 
mode lifetime for all stars. The dashed blue line shows the 1:1 
relation in both panels, and the position of the Sun is marked. 

0.02ppm 2 /iHz _1 . The upper panel of Figure [5] compares 
the proposed relation in Equation (p~2|) . again using as- 
teroseismic masses and radii, as well as effective temper- 
atures from KIC. We observe that the calculated values 
are systematically too high for main-sequence stars, and 
systematically too low for red giants, the la t ter bei ng in 
agreement with the results bv lMosser et al.l (|2011d) . 

In order to compare observed with calculated am- 
plitudes using Equation (fTTj) . estimates of the mode 
lifetimes are required. While it is well established 
that mode lifetimes in giants are significantly l onger 
than in main-sequ e nce stars ([De Ridder et all 120091; 
Chaplin et all 120091; iHekker et al.l l2010al iHuber et al l 
20101: iBaudin et all 120111 ) . t™. is not very well con- 
strained for stars spanning such a large range in evolution 
as considered in our sample. We have therefore decided 
to neglect the influence of r osc in Equation (ITT1) in our 
study. The lower panel of Figure [5] shows the comparison 
between observed and calculated amplitudes using Equa- 
tion (jllj) and assuming solar mode lifetime for all stars. 
We observe that the revised scaling relation works well 
for main-sequence and subgiant stars, but predicts am- 
plitudes that are system atically too h i gh for red giants, 
again in agreement with iMosser et all (|2011d) . 

We note that the upper panel of Figure [5] displays an 
apparent systematic difference in the granulation power 
at Vxaax in the overlapping region of the long-cadence and 
short-cadence sample. This effect is due to a systematic 
difference of the modelled background noise level for stars 
oscillating close to the long-cadence Nyquist frequency 
(283/iHz). We have confirmed that this has a negligible 
effect on the measured amplitudes, as is evident by the 



consistent overlap of datapoints in the bottom panel of 
Figure El 

As noted bv lKieldsen fc Beddind (|2011[ ). the observed 
differences are not unexpected since the relation does not 
take into account potential differences in the intensity 
contrast between dark and bright regions on the stellar 
surface. While our study confirms the relation between 
mode amplitudes and granulation power for stars rang- 
ing from the main-sequence to the red-giant branch, it 
demonstrates that the link between velocity and intensity 
measurements is not yet fully understood for amplitude 
scaling relations. 

4.3. Stellar Activity & Amplitudes 

The offset between observed and calculated amplitudes 
for main-sequence stars noted in Section 4.1 suggests that 
there might be an additional physical dependency which 
is not yet taken into account. Increased stellar activ- 
ity and magnetic fields are known to decrease pulsation 
ampl itudes in the Sun (|Chaplin et alll2000l ; iKomm et al.l 
2000), and have been suggested as a possible cause for 
the u nexpected low amp l itudes of oscilla tions in active 
stars (jMosser et al.l 120091: iDall et al.H201Q[ ). The first di- 
rect evidence for an influence of stellar activity on oscil- 
lation amplitudes for a star oth er than the Sun has been 
found in the F-star HD 49933 (jGarci'a et al.l 12010ft . As 
discussed by KB11, these amplitude changes may be due 
to changes i n the mode lifetimes. 

Recently, IChaplin et al.l (|2011bl ) reported the discov- 
ery for main-sequence and subgiant stars observed with 
Kepler that oscillations are less likely to be detected in 
stars with higher activity. This led them to conclude 
that, just like for the Sun, increased stellar activity sup- 
presses mode amplitudes. Kepler results have further- 
more shown that the detection rate seems to be consid- 
erably lower for subgia n t stars with v max ~ 400—9 00 fiRz 
(|Chaplin et all 1201 lallcl : ISilva Aguirre et all 1201 111 This 
was tentatively explained by increased magneti c activity 
in this stage of stellar evolution, as predicted bv lGillilandl 
( 1985). While it is known that some detection bias due to 
instrumental artefacts mi ght be present in th e frequency 
range from 200-500 //Hz (jGarcfa et ail 120111) . both dis- 
coveries strongly suggest that activity must be consid- 
ered wh en calculating amplitud es. 

While Cha plin et al.l (|2011bf) used a simple detection 
criterion, we attempt here to confirm and quantify their 
discovery by correlating measured pulsation amplitudes 
with stellar activity. Several studies have been devoted 
to stellar activity and variability in the Kepler field 
dBasri et al.l 120101: iDebosscher et all 1201 U iCiardi et all 
120111) with various definitions of activity measures. Since 
our study covers a large range of oscillation timescales 
and amplitudes, we need to define an activity measure 
which separates variability due to oscillations from vari- 
ations due to activity, and at the same time remains sen- 
sitive to local minima and maxima in the light curve. 
To do this, we first smoothed each light-curve with a 
quadratic Savitzky-Golay filter (Sav itzkv fe Golavlll964|) 
with a full-width corresponding to 20 times the deter- 
mined oscillation period, and then measured the abso- 
lute maximum deviation of the smoothed curve from 
its mean. We denote this measure of activity as r sg . 
Note that we have repeated the analysis with smoothing 
widths ranging between 10-50 times the oscillation pe- 
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Fig. 6. — Amplitude versus f ma x showing only stars observed 
in short cadence, with the logarithm of the activity range color 
coded. Red and blue points correspond to high and low activity, 
respectively. The position of the Sun with a color corresponding 
to the mean solar activity is also shown. 



riod, but found not significant difference in the results. 
When evaluating r sg , it must be kept in mind that our 
sample includes light curves with two very different time- 
bases (^30d for short cadence, and ^500 d for long ca- 
dence). The samples therefore probe different activity 
time scales, and have hence been separated in the subse- 
quent analysis. 

To correct our activity measure for the influence of 
shot noise, we performed the following simulations. For 
each target light curve, we produced synthetic time- 
series including a sinusoidal variation with an ampli- 
tude corresponding to the measured r sg value. We then 
added white noise with a standard deviation correspond- 
ing to the apparent magnitude of the target, which 
was estimated from th e mi nimal noise le v els gi ven in 
IGilliland et all (|2010ah and I Jenkins et all (12010) . We 
then measured r sg for the synthetic light curve as de- 
scribed above, and calculated the difference between this 
value and the value of r sg measured from the noise-free 
time series. This process was repeated 500 times for each 
star, and the median of the resulting distribution was 
taken as the correction value for the observed r s<r . In 
summary, it was found that the typical correction val- 
ues were about 8% for short-cadence data, and < 1% for 
long-cadence data. 

Figure |5] shows the amplitude-f max relation for the 
short cadence sample, with the logarithm of the activity 
measure color-coded. We see that the most active stars 
(shown in red) have generally lower ampli tudes than less 
active s tars (shown in blue). As found bv iChaplin et al.l 
(|2011bf ). the Kepler sample is generally less active than 
the average Sun (r sg: = 0.7ppt), in agreement with 
the bias towards lower predicted amplitudes found in the 
previous section. To quantify this, we first removed the 
dependence of oscillation amplitudes on stellar proper- 
ties by dividing them by the calculated amplitudes using 
Equation (J9j) with the best-fitting coefficients derived in 
the previous section. The resulting normalized ampli- 
tudes are plotted versus the logarithm of the activity 
measure in Figure [7] 

Figure [7| is split into three ranges of v msK , roughly di- 
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Fig. 7. — Upper panel: Residuals of observed amplitude divided 
by the calculated amplitudes using Equation[9]versus the logarithm 
of the activity measure for all stars with v max > 1000/^Hz. The 
red solid line shows an unweighted linear fit. Middle panel: Same 
as the upper panel but for all stars with f max = 300 — 1000 /xHz. 
Lower panel: Same as upper panel but for all long-cadence stars 
Kai < 300/xHz). 



viding stars on the main-sequence (top panel) , subgiants 
(middle panel) and red giants (bottom panel). The influ- 
ence of activity on amplitudes appears the strongest for 
subgiant stars, followed by a weak correlation for main- 
sequence stars and no visible correlation for red giants. 
Unweighted linear fits to the linear-log plots yield slopes 
of -0.13±0.04, -0.24±0.07 and 0.00±0.02, respectively. 
While the correlations for unevolved stars are formally 
significant, a thorough comparison with other methods 
has shown that the correlation for subgiant stars is only 
confirmed in three out of six methods, while the results 
for main-sequence and red-giant stars are consistent with 
zero within 3c for all methods. Although this compar- 
ison cautions us not to draw any definite conclusions, 
we note that the exponential decrease of amplitude with 
increasing stellar activity for subgiants would be in-line 
with both observational and theoretical evidence for en- 
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log (T cH /K) 

Fig. 8. — %ax versus effective temperature for all stars in our 
sample. Color codes refer to the logarithm of the pulsation am- 
plitude, with blue and red colors marking the highest and lowest 
amplitudes, respectively. The dashed line shows the empirically 
determined cool edge of the instability strip from observations of 
8 Scuti stars. The theoretical cool edge for the fundamental radial 
mode of 8 Scuti pulsations is shown for a 1.7 Mm model bylfloudekl 
(2000) (solid circle) and a range of masses by Duprct ct al. (2005) 
(dashed-dotted line). Black solid lines are solar- metallicity BaSTI 
evolutionary tracks with masses as indicated in the plot. The star 
symbol indicat es the first detection of hybrid solar-like and 8 Scuti 
oscillations by Antoci ct al. (20fJ). 

hanced magnetic act ivity in these stars (iGillilandl 119851 : 
Chapl in" et alll20~llaf) . A confirmation of the results pre- 
sented here will have to await longer timeseries (in par- 
ticular for short cadence data) which will allow a better 
estimate of the activity measure and reduce the uncer- 
tainty in the adopted stellar properties through better 
constraints on ^ max and Aza This will then allow a more 
in-depth investigation of additional physical influences on 
oscilla tion amplitudes suc h as stellar activity and metal- 
licity (|Samadi et al.ll2010fl . 

5. THE COOL EDGE OF THE INSTABILITY STRIP 

The cool edge of the instability strip is widely be- 
lieved to be the dividing line between coherent pulsations 
driven by the opacity (k) mechanism and solar-like os- 
cillations driven by convection. Considerable work has 
been devoted to establishing this boundary both empiri- 
cally using S Scuti and 7 Doradus stars (see, e.g.. iBregerl 
I1979L iPamvatn vkh 200(J as well as the oretically (see, 
e.g.. Houdek et alj|1999bt IHoudeM l2000t IXiong fc Deng] 
120011; IDupret et all 120041 |2005J). The large number of 
stars for which Kepler has detected oscillations now al- 
lows the first test of this boundary from "the cool side" , 
using solar-like oscillations. 

Figure [8] shows all stars of our sample in a modi- 
fied H-R diagram, in which we have replaced luminosity 



with 1/^max- We also show the empiricall y determined 
cool e dge of the instability strip taken from IPamvatnvkhl 
(2000) as a dashed line, as well as the theoretical cool 
edge for the fundamental radial mode of S Scu t i pul - 
sations for a 1.7 Mq model taken from iHpudeld (l2000f) 
(filled circle) and for a range of masses by IDupret et all 
(2005) (dashed-dotted line). To convert the lines from 
the log L — log T e g to the log v max — log T c s plane, we de- 
termined the closest matching grid point of evolutionary 
tracks with different masses to the log L — log T e g rela- 
tion, and then fitted a straight line to the resulting points 
in the log f max — log T c g plane. T his has been done both 
for B aSTI and a set of ASTEC (|Christensen-Dalsgaardl 
2008) evolutionary models, and both sets of models 
yielded consistent resu lts. The emp i rically determined 
cool edges taken from IPamvatnvkhl (|2000n in Figure [8] 
are: 

log(T eS /K) = -0.045 log(L/L ) + 3.893 (13) 



log(T eff /K) = 0.048 log(i/ max //xHz) + 3.702 . (14) 

Equation (|14p appears in good agreement with the 
hottest stars for which solar-like oscillations are observed. 
However, it is important to note that the recent work by 
iPinsonneault et all (|2011l ) and iMolcnda-Zakowic z et all 
(|2011D suggest that many stars are significantly hotter 
than the effective temperatures given in KIC, with dif- 
ferences of up to 250 K. With this in mind, the fact that 
we observe stars close to the theoretical limit is evidence 
that the separation between the classical and stochastic 
pulsations is probably not a sharp dividing boundary. 
Indeed, we indeed do not observe a gradual decrease in 
amplitude towards the cool edge of the instability strip 
and the recent discove ry of solar-lik e and 5 Scuti oscilla- 
tions in HD 187547 bv lAntoci et~aTI (|2011l) using Kepler 
data has confirmed the existence of hybrid pulsators (see 
star symbol in Figure [8]). Interestingly, the observed 
amplitude for this star is roughly a factor four higher 
than expected from the traditional scaling relations dis- 
cussed in this paper. This was tentatively explained by 
the increased mode lifetimes which are roughly 4-5 times 
higher than expected from scaling relations based on ef- 
fectiv e temperature (jChaplin et al.l I2009L iBaudin et all 

unni). 

To test the cool edge of the instability strip more quan- 
titatively, we must account for selection effects. Figure 
[5] shows a close-up of the main-sequence region of the 
HRD, this time plotting luminosity on the ordinate and 
including stars for which we did not find evidence for os- 
cillations. Luminosities for stars without detections have 
been calculated based on values in the KIC assuming a n 
uncertainty on the radius of 40% (jVerner et al.ll2011bj) . 
while for stars with detections Equations (TTJ) and ^ were 
used. Although the error bars are large, we observe an 
excess of hot stars for which no oscillations have been 
detected, indicating that the observed cut-off in detec- 
tions may not entirely due be due to selection bias. We 
therefore conclude that the empirical red-edge given in 
Equations (ITU)) and (fH| gives a good approximation for 
the transition between opacity driven and solar-like oscil- 
lations. Further work will be needed to quantify this ob- 
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3,88 3.86 3-84 3.82 3,80 3,78 3.76 
log (T eH /K) 

Fig. 9. — H-R diagram of the low-luminosity region shown in 
Figure [Sj comparing stars for which no oscillations have been de- 
tected (open black diamonds) with our sample of detections (red 
filled diamonds). A typical error bar for both samples is shown in 
the top region of the plot. The dashed line, dashed-dotted line and 
filled circle are the same as in Figure [8] The star symbol indicates 
the first detection of hybrid solar-like and 8 Scuti oscillations by 
lAntoci et all <20Tll) . 

servation and compare the results with theoretical mod- 
els of excitation and damping of solar-like oscillations, 
as well as possible further detections of hybrid solar-like 
oscillations and classical pulsation in hot stars. 

6. CONCLUSIONS 

We have studied global oscillation properties in ^1700 
stars observed by Kepler to test asteroseismic scaling re- 
lations. Our main findings can be summarized as follows: 

• By comparing evolutionary models with observa- 
tions, we have shown that the scaling relations for 
i/ max and Av are in qualitative agreement with ob- 
servations for evolutionary stages spanning from 
the main-sequence to the He-core burning phase 
of red giants. The difference in the Ai/-^ max re- 
lation between evolved and unevolved stars can be 
explained by different distributions of effective tem- 
perature and stellar mass, in agreement with what 
is expected from the scaling relations. A more 
quantitative test of scaling relations for v max and 
Av will have to await the determination of funda- 
mental properties from independent methods such 
as spectroscopy and long-baseline interferometry. 
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• We have shown that (L/M) s scaling for oscillation 
amplitudes fails to reproduce the amplitude-zy max 
relation for red giants. We have verified that a 
revised scaling relation using a separate mass and 
luminosity dependence reproduces the observations 
better, and a relation with L s and M* coefficients of 
s = 0.838±0.002 and t = 1.32±0.02 matches ampli- 
tudes for field stars ranging from the main-sequence 
to the red clump to a precision of 25% when adopt- 
ing stellar properties derived from scaling relations. 
The calculated amplitudes for main-sequence and 
subgiant stars in our sample are systematically un- 
derestimated by up to 15%, indicating that there 
might be an additional physical dependence which 
is not yet taken into account. 

• We have investigated the connection of stellar ac- 
tivity with the suppression of oscillation ampli- 
tudes in main-sequence, subgiant and red-giant 
stars. We find evidence that the effect is strongest 
for subgiant stars, but caution that these results 
will have to await confirmation with longer time 
series providing a better estimate of the activity 
measure and reduced uncertainties on stellar prop- 
erties used to calculated amplitudes. The present 
data do not yield strong evidence that stellar activ- 
ity contributes significantly to the underestimation 
of calculated amplitudes for main-sequence stars. 

• We have investigated the cool edge of the instability 
strip using detections of solar-like oscillations. We 
find good agreement with the empirically and the- 
oretically determined cool edge using 8 Scuti stars, 
but note that many stars showing solar-like oscil- 
lations may overlap with cooler 5 Scuti stars, in 
agreement with the recent first discov ery of hybrid 
solar-l ike and 5 Scuti oscillations by lAntoci et al.l 
(l20ilr) . 
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